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The results of a recent quenched lattice QCD simulation for charmed baryons are presented. In contrast 
to a previous calculation, it is found that hyperfine splittings are in agreement with quark model expectations 
and comparable to experimental values. Preliminary calculations using lattice NRQCD yield results which are 
compatible with those obtained using a Dirac- Wilson action of the D234 type for the charm quark. 



There have been many applications of lattice 
QCD to heavy-flavoured hadrons but most of the 
work has been restricted to the meson sector. The 
most complete study of heavy baryons in a rela- 
tivistic framework is the one done in the UKQCD 
collaboration^ using an O(a) improved fermion 
action in quenched approximation. In that work 
the overall features of the baryon spectrum were 
found to be in agreement with experiment but 
hyperfine splittings in both charm and bottom 
baryons were found to be considerably smaller 
than those predicted by phenomenological mod- 
els M and observed experimentally. In fact, the 
central values for hyperfine splittings found in M 
were negative. Such a result is difficult to under- 
stand and would pose a severe problem for lattice 
QCD. Therefore further investigation seems jus- 
tified. 

Due to limitations in computing resources it 
was not possible in this work to use the same lat- 
tice spacing and volume as used in the UKQCD 
calculation. Rather, a special strategy had to 
be adopted. We work on a more coarse lattice 
(^0.2fm) with a highly improved action. Past 
experience has shown that results of reasonable 
accuracy may be obtained with such lattices |3_). 
In order to check the calculation, the spectrum of 
baryons in the light quark (u,d,s) sector was cal- 
culated at the same time. As well, meson masses 
for both heavy and light quarks were calculated. 
The results of all these calculations (see jl for 
details) are in reasonable agreement with exper- 
imental values and with the results obtained at 
small lattice spacing flM. 



The explicit form of the gluon and fermion ac- 
tion may be found in [|l| . The fermion action is of 
the D234 type|7| with tadpole improvement im- 
plemented using the Landau link. The lattice was 
anisotropic with a bare aspect ratio a s /at of 2. A 
total of 420 quenched configurations on a 10 3 x 30 
lattice were analyzed. 

The results for singly charmed baryons are 
given in Table 0. In contrast to M the hyper- 
fine splittings are positive and large in magnitude. 
Figure p] shows the mass difference S* — E c as 
well as the results for the E-hyperon and nucleon 
- delta compared to experimental values. Also 
shown are the results of lattice calculations [EM 
done with different actions on lattices with much 
smaller lattice spacing than used here. The re- 
sults of our simulation are seen to been quite com- 
patible. 

The masses of doubly charmed baryons 
have also been calculated. The results are 
3.598(13)GeV and 3.682(20)GeV for S cc and E* c 
respectively. For Q cc and f2* c we get masses of 
3.697(10)GeV and 3.775(12)GeV. The hyperfine 
splittings are only slightly smaller than for singly 
charmed baryons. 

In lattice QCD masses are extracted from cor- 
relation functions of hadron operators. One dif- 
ference between the present calculation and that 
of Rcf.|l| is in the choice of baryon operators. In 
fjj calculations are done with interpolating oper- 
ators of the form: 
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where 1,1' are light quark fields and h is a heavy 
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Figure 1 . Spin splitting for baryons with different 
flavour versus inverse of the corresponding pseu- 
doscalar meson mass. The open squares are ex- 
perimental values and the triangles are the re- 
sults of this work. For comparison the results of 
|5) (filled squares) and ||] (filled circles) are also 
shown. 



Table 1 

Masses of singly charmed baryons. Masses are 
given in GeV, mass differences are in MeV. The 
experimental values are taken from g] except for 
S' which is from |B| 





This work 


UKQCD[1] 


Experiment 


A c 


2.304(39) 


2.27 It +f 


2.285(1) 


s c 


2.465(22) 


2.46 ±1 ±1 


2.453(1) 


E* 


2.557(30) 


2.44 +1 +4 5 


2.518(2) 
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2.454(21) 


2.41 +3 +4 


2.468(2) 


'— 'c 


2.579(14) 


2.51 ±1 !« 


2.575(3) 


1— 'c 


2.672(16) 


2.55 ±| ±b 


2.645(2) 


n c 


2.664(12) 


2.68 ±5 i§ 


2.704(4) 


n* c 


2.757(14) 


2.66 ±1 t 6 7 




£*-E c 


91(25) 


17 +12 +3 
11 -31 -2 


65(2) 


1— 'c i ~ 'c 


94(13) 


on +12 +2 
~ zu -24 -3 


70(4) 


ft* -Q c 


94(10) 


90 +6 +3 
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field (for M , I' may be the same flavour as I). 
This choice of operators is motivated by heavy 
quark symmetry where baryons are classified ac- 
cording to the total spin (0 or 1) of the light 
flavours. 

Our approach is different. The aim is to do a 
unified analysis, incorporating all flavours (quark 
masses) in a common framework. We therefore 
start with the nucleoli operator for which a com- 
mon choice ]10| , [il| is 



'(ulCj 5 d b )u c 



(2) 



where u and d are up and down quark fields. It 
can be shown that this operator contains the com- 
ponent 



abc 



(wI c '7m'W&)757 m 4 



(3) 



which is expected to survive in the limit where 
the mass of the unlike quark (d in this case) is 
taken to be very heavy. Therefore the familiar op- 
erators used in light baryon spectroscopy should 
also be applicable, with appropriate substitution 
of fields, in the heavy baryon sector. The forms 
used may be found in jl| . 

In principle, masses should be independent of 
the interpolating operators used as long as the 
quantum numbers are the same. In view of the 
difference between the results of the present cal- 
culation and Ref.M it is worthwhile to check if 
this independence of operator is obtained in prac- 
tice. We have therefore analyzed a subset of our 
configurations using operators of the form (1) and 
compared the results with those obtained using 
the operators given in (3. Typical comparisons 
are shown in figure 0. The quantity shown in the 
figure is M(t) = ln(G(t)/G(t + 1)) where G(t) 
is the hadron correlator and for large times t it 
is equal to the mass. It is seen that the mass is 
indeed the same for the two different choices of 
interpolating operator. 

In the continuum limit the D234 action de- 
scribes Dirac fermions. For heavy quarks on the 
lattice, it can be advantageous to use a nonrela- 
tivistic approach |12j|. For b-quarks this seems to 
work well but for charmed quarks it is less clear 
whether lattice NRQCD is adequate US A 
particular concern is that hyperfine splittings in 
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Figure 3. The mass difference AM for f2* — il c 
calculated using lattice NRQCD for the charmed 
quark. The intersection of the cross hatched 
bands indicates the result obtained with the D234 
action. 
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Figure 2. M(t) calculated using two different 
baryon operators. The upper figure is for f2 c and 
the lower is for SI*. The results using the op- 
erators (1) are shown by circles, those using the 
operators of [0 are shown by triangles. 



the meson sector seem to be underestimated com- 
pared to experimental values |13|-[l5[l . Therefore it 
is of interest to investigate the baryon sector. 

We have carried out a preliminary analysis of 
charmed baryons using lattice NRQCD, instead 
of the D234 action, for the charmed quark. The 
result for the mass difference Q* — fi c is plotted 
in Figure and for fl* c — fl cc in Figure 0. For 
comparison, the results obtained using the D234 
action for the charmed quark are also shown. Cal- 
culations were done for two different values of the 
nonrelativistic charm quark mass which should 
be tuned to reproduce the charmed baryon mass 
(indicated by the vertical cross hatched area in 
the figures). Interpolating between the two points 
calculated using NRQCD, one sees that the result 
for the spin splitting using NRQCD is compatible 
with that using the D234 ("relativistic") action. 

To summarize, we have calculated the spec- 
trum of singly and doubly charmed baryons in 
quenched lattice QCD. In contrast to a previous 
calculation Q, it was found that hyperfine split- 
tings are positive as expected, for example, from 
the quark model and are comparable in magni- 
tude to experimental values. This result is quite 
robust. It is independent of the choice of inter- 
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Figure 4. The mass difference AM for Sl* c — fi cc 
calculated using lattice NRQCD for the charmed 
quark. The intersection of the cross hatched 
bands indicates the result obtained with the D234 
action. 



polating operators for the baryons. Preliminary 
results also indicate that use of an NRQCD de- 
scription for the charmed quark yields values for 
the spin splittings compatible with those obtained 
using the D234 action. 
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